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^ . Abstract 

Q^l The method of Mellin-Barnes representation is used to calculate di- 

mensionally regularized massive on-shell double box Feynman diagrams 
contributing to Bhabha scattering at two loops. 



1 Introduction 



The method of MeUin-Barnes representation has turned out to be very successful 
for the analytical evaluation of two-loop Feynman diagrams with four external lines 
within dimensional regularization pP . The problem of the evaluation of massless on- 
shell double boxes was solved in 013111 El, with multiple subsequent applications 
to the calculation of two- loop scattering amplitudes in gauge theories [B] . First cal- 
culations of massless on-shell four-point three-loop Feynman integrals were done in 

m 

Complete algorithms for the evaluation of massless on-shell double boxes with 
one leg off shell were also constructed. Master integrals were calculated using the 
MB representation (first results in 0) and the method of differential equations (DE) 
[Oj, where a systematic evaluation was described in [10]. The reduction to master 
integrals was done using a Laporta's idea jTTj in ^0]- AH results are expressed in 
terms of two-dimensional harmonic polylogarithms jTO] which generalize harmonic 
polylogarithms (HPL) ^21- Let us stress that this very combination of reduction 
based on [TTl [TU] and DE was successfully apphed in numerous calculations, e.g. to 
various classes of vertex diagrams ^3J ^^ . 

Another important class of Feynman integrals with one more parameter, with 
respect to the massless on-shell double boxes, are massive on-shell double boxes 
relevant to Bhabha scattering at NNLO. Bhabha scattering plays an important role 
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Figure 1: Massive on-shell double boxes: (a) planar double box of the first type, 
(b) planar double box of the second type and (c) non-planar double box. The solid 
lines denote massive, the dotted lines massless particles. 



as a luminosity monitor for e+e~ colliders. Because of this phenomenological interest, 
the two-loop corrections to Bhabha scattering have been studied extensively in recent 
years. In ^l^, the complete two- loop virtual QED corrections were calculated in 
the massless fermion approximation. The infrared singularity structure of this result 
was studied in [TBI. The first result for a two-loop master double box with nonzero 
electron mass was obtained in J7]. Planar two-loop box diagrams with one-loop 
insertion and without neglecting the fermion masses were calculated by DE in |18j . 
This work was continued to calculate the subset of all two-loop diagrams involving a 
closed fermion loop in J^. The two-loop contribution which factorizes into squares 



of one- loop diagrams was calculated in [201 ■ Numerical results for the three master 
integrals shown in Fig.^ were given in J2I] for Euclidean points. Despite all these 
efforts, a complete two- loop differential cross section is not available yet [22j, mainly 
because the most complicated master integrals, the two-loop box diagrams, have not 
all been calculated so far. 

2 Evaluating by Mellin-Barnes representation 

2.1 Planar double boxes of the first type 

The general double box Feynman integral of the first type (see Fig. la) takes the form 

BpL^i{ai,...,a8;s,t,m ; e) 



(P - m2)«i [{k + pi)2)]"2 [(k + pi+ P2Y - m'^Y'' 



1) 



[(/ + P1+P2Y - m?)Y^ [{I +P1+P2+ PzYY" {P - m?Y^ [{k - /)2]"7 ' 

where s = {pi +^2)^? t = {p2 +^3)^, and k and / are the loop momenta corresponding 
to the left and the right box. Usual prescriptions k"^ = k'^ + iO, s = s + iO, etc. are 
implied. For convenience, we consider the factor with {k+pi+p2+P3Y corresponding 
to the irreducible numerator as an extra propagator but, really, we are interested only 
in the non-positive integer values of ag. 

A first analytical result for the planar double box of the first type, was obtained 
in ^\ using the following sixfold MB representation: 

BpL,iiai,...,a8;s,t,m'^;e) - 



n,=2,4,5,6,7r(a,)r(4 - a4567 " 2e){-s)-''+ 



-4+2e 



5 /^2\^l+^5 



1 /■+'°°j A^ m^Y'^'' ftyT{a2 + w)T{-w)T{z2 + z^)T{z3 + Z4) 



X- —I dwYldzj. , . , 

{271%)'^ J~ioo jj-^ \~^J ^^^ T[ai + Z3 + Z4)T[a3 + Z2 + Z4) 

r(4 - ai3 - 2a28 - 2e + ^2 + 23)r(ai238 - 2 + e + z^ + Z5)T{a7 + w - Z4,) 

r(4 — 046 — 2a57 — 2e — 2w — 22:1 — Z2 — z^) 
r(a.4567 -2 + e + w + Zi- -24)r(a8 - ^2 - ^3 - Z4)T{-w - 2:2 - -23 - Z4) 



X 



r(4 - ai238 -2e + w- Z4)T{as - w - Z2 - z^ - z^) 
r(a5 + w + Z2 + z^ + z^)V{2 - a^Qr - e - w - zi - Z2) 
r(4 - an - 2a28 - 2e + Z2 + Z3 - 2z^) 
xr(2 - a457 -e-w- zi~ Z'i)T{2 - ai2s -t + Z2- z^)T{2 - 0238 -t + z^-zr^) 
xr(4 - a46 - 2a57 -2e-2w-Z2- Z3)V{-z{)V{-z^) , (2) 

where a = ai...8, 04567 = 04 + ^5 + 06 + O'T: ^Hs = ^i + %, etc., and integration contours 
are chosen in the standard way. 



To evaluate the master double box BpL^i{l, ... ,1,0; s,t, w?; e) the standard pro- 
cedure of taking residues and shifting contours [2^ was used [TTj, with the goal to 
obtain a sum of integrals where one may expand the integrands in a Laurent series in 
e = (4 — (i)/2 (where d is the space-time dimension within dimensional regularization 
PP). Then one can use the first and the second Barnes lemmas and their corollaries to 
perform some of the MB integrations explicitly. In the last integrations which usually 
carry a dependence on the external variables, one closes the contour in the complex 
plane and sums up the corresponding series. (See j21j for details of the method.) 

In |22] it was announced that the massive on-shell double boxes with two re- 
duced lines have been calculated by DE, with a solution of reduction problem by the 
method of ^TJOni- However an extension of these results to all master integrals meets 
problems. Differential equations of third order and higher are encountered there. 

On the other hand, the method of MB representation can certainly be applied 
to arbitrary massive on-shell double boxes. To illustrate this point, we present the 
evaluation of typical complicated master integrals from this class. 

As has been demonstrated in j3], it can be useful to use a double box with a 
numerator as a master integral in order to avoid coefficients of order 1/e in the course 
of the tensor reduction when calculating a scattering amplitude. Therefore we present 
here results for the planar double box BpL^i{l, . . . , 1, —1; s, t, rn?; e). 

Our calculation of this Feynman integral is quite similar to the previous case 
[T7] . It is based on (J2I) . After a preliminary analysis of the integrand one observes 
that it is reasonable to start the procedure of resolution of the singularities in e, 
by shifting contours and taking residues, from taking care of the gamma function 
r(— 1 — Z2 — z^ — Z4). Here one meets the same problem as in the massless case (see 
the third reference in j3]) connected with spurious singularities in MB integrals. It 
can also be cured in the same way, by introducing an auxiliary analytic regularization. 
To do this, we have chosen ag = —1 + x. Then the singularities in the MB integrals 
are first resolved with respect to x and then with respect to e when x and e tend to 
zero. With this complication, the key gamma function becomes r{x — 1 — Z2 — z^ — Z4) . 
As a result of the procedure of resolution of the singularities in x and e, the spurious 
singularities in x drop out, and we are left only with six MB integrals which can be 
evaluated by expanding the integrands in e and evaluating the corresponding finite 
MB integrals. 

We have obtained the following result: 



5pL,i(l,...,l,-l;s,t,m ;e 



s^fm^)^*^ 



^ + ^ + 4i(x) + 4(x) +41nl±^4(x,y) + 0(6) 



with 



4W - I'-'l^. (3) 



c| X = L3 ^ ^ - 2n + 2n2-n^r L2 7 br 

+2 In \—^ Lis f ^^l + ^ ln=^(l - x) + ^ \n\l - x) ln(l + x) 
1+x V2/0 2 

— In(l-x)ln2(l + x) + -ln^(l + x) -2 In 2 ln(l-x)ln^— ^ 
2 6 1 + X 

-yl^5-21n^[^ln^-y(31n(l-x)-21nx-ln(l + x)) 
+ ln221ni^ + ^ln2-C(3), (4) 

4(x) = 8i7(-l,0,0,l;^^^) , (5) 



where x = l/yl — Am? / s and H(—l, 0, 0, 1; z) is a HPL 

The two other contributions to the finite part are more cumbersome: 

0^2(3;) = 3Li4((l-x)V(l + a;)^) -8Li4((l-x)/2)-8Li4(-2x/(l-x)) 

-12Li4 (-(1 - x)/(l + x)) - 4Li4 ((1 - x)/(l + x)) - 8Li4 (2x/(l + x)) 

-8Li4 ((1 + x)/2) - 2^2,2 ((1 - x)V(l + xf) - 2(4/2 - m, - 3p,)Li3 ((1 - x)/2) 

+4(m^. - p^)Li3 (1 - x) - 2(2/2 + h - m^ - p^)U3 (-4x/(l - x)^) 

+4(m^. - pjLis (-x/(l - x)) + 2(2/2 + 2/^ - 3m^. + p^)Li3 (-2x/(l - x)) 

-2(3/2 + /. - m, - 2p,)Li3 ((1 - xf/{l + xf) 

-2(2/2 + lx- m^ - pjLia (4x/(l + x)^) - 4(m^ - p^)Li3 (1/(1 + x)) 

-2(6/2 + 2/, + m, - TpjLis ((1 - x)/(l + x)) 

-4(5/2 + 3/, - 2m, - 3p,)Li3 (-(1 - x)/(l + x)) - 4(m, - p,)Li3 (x/(l + x)) 

+2(2/2 + 2/, + m,- 3p,)Li3 (2x/(l + x)) - 2(4/2 - 3m, - p,)Li3 ((1 + x)/2) 

+4(/2 + /, - m,)Li3 (-(1 + x)/(l - x)) + 4Li2 ((1 - x)/2)' + 2Li2 (-x)' 

+ (m, -Px){-Ql2 - 12/, + 3m, + 3p, + 21n(-s/m^))(Li2 (-x) - Li2 (x)) 

-4Li2 (-x) Li2 (x) + (2/2 + /, - 2m,) (m, - p,)Li2 (-4x/(l - x)^) 

H^x - Px)(2/2 + 4/, - m, - ln(-s/m2))Li2 (-2x/(l - x)) + (2/^ + 2/2/, 

+II/2 - 4/2m, - 3/,m, - m^ + /,p, + 6m,p, - 3p2)Li2 ((1 - x) V(l + xf) 

-2(2/2 + k- 2p,){m, - Px)U2 (4x/(l + x)=') - (4/2 + 4/2/, + /^ - 8/2m, 

-4/,m, -ml- 7rV3 + 10m,p, - 5p^)Li2 ((1 - x)/(l + x)) - {Alj + Akk + ll 

-Ahm^ - 2ml - tt^ - 4/2P, - 4/,p, + 8m,p, - 2pl + 4Li2 ((1 - x)/(l + x)) 

-2Li2 (-(1 - x)/(l + x)))Li2 (-(1 - x)/(l + x)) 

-(m, -p,)(2/2 + 4/, -m, - 2p^ - ln(-s/m2))Li2 (2x/(l + x)) 

+6Li2 (-x) Li2 ((1 + x)/2) + 2(-3m, - 4/,m, + 2m^. + 3p, + 4/,p, - 4m,p, + 2pl 



+2m^ ln(-s/m2) - 2p^ \n{-s/m^) - SLia (a;) + 2Li2 ((1 + x)/2))U2 ((1 + x)/2) 

-2(3m^ - 8/2"^^ - lO/^m^. + 2m^ - 3p^ + 8/2^^ + lO/^Px + 4m^p^ - 6p^ 

+2m2: ln(— s/m^) — 2px ln(— s/m^) + 3Li2 {—x) — 3Li2 (x) 

+4Li2 ((1 + a;)/2))Li2 ((1 - x)/2) + 2Li2 {xf 

+2/2 - Qllm^ + 4/2^^ + Qlll^m^ + 6/2m^ - 41/2^^ 

-32/2/xm^ + 2/^m^ + (65/2m;^)/3 + (23/^m^)/3 - (49m^)/8 - (2/^7r2)/3 

+m^7r2 + (5/2m^7r2)/3 - {l^m^^'^)/2 - {7ml^'^)/l2 + (377r^)/360 

+6/2PX - I2/2PX - 6/2/xPx + 70llm^p^ + 52l2lxm^Px 

-AllrrixPx - Qmlp^ + Ukmlp^ + ISl^mlp^ + {4mlpx)/3 

-ttV - (^27rV)/3 + (/^7rVx)/2 - {llm.^7r'^p^)/6 - U2PI 

-1711pI - 20/2/xPx + 2/^P^ + Qm^pl - 87l2m.^pl - 43/^m^p^ 

-{Almlpl)/A + i7n'pl)/A + (UlkpD/S + i67lxpl)/3 + in3m.,pl)/3 

-{165pl)/8 - 2(2/2 + 4/, -m^- p,){m^ - p.f \n{-s/m^) 

+2{m^ - p^f \n^{-s/m^) + 2(m^ - p^f \n^{-t/m^) 

+I2/2C3 + 6/.C3 - 4m,C3 - 8P.C3 , (6) 

4(x,t/) = 2(Li3((l-x)/2)-Li3((l + x)/2) + Li3((l-x)y/(l-xy)) 

-Lia (-(1 + x)yl{l - xy)) + Lig (-(1 - a;)y/(l + xy)) - Lig ((1 + x)yl{l + a;t/))) 

+Li3 ((1 + x){l- y)/(2(l - xy))) - Li3 ((1 -x){l+ t/)/(2(l - xy))) 

-Li3 ((1 - x)(l - y)/(2(l + xy))) + Li3 ((1 + x)(l + t/)/(2(l + xy))) 

+2(/2 - m,)Li2 ((1 - x)/2) + 2(p, - /2)Li2 ((1 + x)/2) - (m,, - p,,)Li2 ((1 - y)/2) 

+2(m^j^ - Pxy){U2 i-y) - Li2 (y)) + {m^y - Py)U2 (((1 + x)y)/{l + y)) 

+ {m^y - Pxy)'^k ((1 + y)/2) + (m^y - Py)Li2 ((1 - xy)/{l + y)) 

-{Pxy - Py)'^i2 {{y - xy)/{l + y)) - 2{ly + m^- m^y)U2 ((1 - x)y/{l - xy)) 

+2{ly - m^y + Px)U2 (-(1 + x)y/{l - xy)) 

-2{ly + m^ - Pxy)U2 (-(1 - x)y/{l + xy)) + 2{ly + Px - Pxy)^h ((1 + x)y/{l + xy)) 

-{Pxy - Py)^i2 ((1 + xy)/{l + y)) 

+ {h + m^y -my- Px)lji2 ((1 + x)(l - y)/{2{l - xy))) 

-{h -rn.^ + m^y - Py)'^k ((1 - x)(l + y)/(2(l - xy))) 

-{h -m^ ~ my + p^y)U2 ((1 -x)(l - y)/(2(l + xy))) 

+ (/2 - Px + Pxy - Py)U2 ((1 + x)(l + y)/(2(l + xy))) 

-{m^y -my+ p^y - Py){2U2 (x) - 2Li2 (-x) + U2 (-2x/(l - x)) 

-Li2 (2x/(l + x))) + {2m^m^y - 2m^my - m^yrriy - 2m^yp^ + 2myp^ 

-f-ZTJlxPxy ~r miyPxy Z,PxPxy ^^'xPy ~r mL^yPy + Z,PxPy PxyPy )''2 ~r '^''yTTLxTTLxy 

+lymly - mainly - 2lymxmy + mxirixymy - mlyiriy + (m^j^m^)/2 



-2lym^yP^ + 2mlyP^ + 2lymyp^ + m^yniyPx 

+ (57r2(m^j, - Pxy))/Q + ^lyiTi^Pxy - m^myPxy - {mlpxy)/2 - 2lyPxPxy 

-myPr^Pxy - lyPly " 2771^^^^^ + TTlyply + P^ply " 2lymxPy " 2lymxyPy " m^m^yPy 

-tZTflxTlilyPy -\- ZlyPxPy ^f^xyPxPi/ ^T^yPxPy ~\~ ZlyPxyPy + 'jf^xPxyPy ~\~ PxPxyPy 

-rrixpl + {3mxypl)/2 + pxpl - {?>Pxypl)/2 , (7) 



where y = l/y 1 — -iiv? /t and 5*2,2 is a generalized polylogarithm |2S]. The following 
abbreviations are also used: Cs = C(3), /^ = Inz for z = x,y,2, pz = ln(l + z) and 
TJiz = ln(l — z) for z = x,y, xy. 

The results have been checked numerically in two independent ways: The finite 
MB integrals, obtained after the resolution of singularities, have been evaluated by 
numerical integration along the imaginary axis. The parameter space for the Monte 
Carlo integration was maximally 4-dimensional in this case. In addition, the overall 
coefficients of the poles have been calculated at Euclidean points by the completely 
independent method of sector decomposition [2Sl I2I]- 

2.2 Planar double boxes of the second type 

For the general planar double box of the second type (see Fig. lb) 

d'^A;d'^/ 



BpL,2iai,...,a8;s,t,m ; e) 



(P - m2)«i [(A; + pi)2)]«2 [{k + pi + P2Y - m^J^s 
[{k + pi+p2+P:if]~^^ ,„, 



[(/ + pi+ P2Y)Y^ \{l +Pi+P2 + P-iY - m?Y^ (/2)«6 [(A; - /)2 - rn?Y^ 
a sixfold MB representation can be derived j^H] similarly to the first case: 

5pL,2(ai,...,a8;s,t,m^;e) 

m'^lA^ i-lY 



n,=2,4,5,6,7r(a,)r(4 - a4567 " 2e) (-s)"-4+2e 

1 /■+^°° A . /^!^V''''Vi V A rr- ^ ^(a2 + ^l)^(a4 + ^2 + -g4) 

''{27^ifJ-.oo}},'\-s) Vs) }X^ ''^ T{a,-z,)T{a,-z,) 

r(4 - 0445667 - 2e - 2:2 - ^3 - 2zi)T{aQ + Z3 + Zj) 

r(4 - a445667 -2e- Z2- z-i- 2zi - 22;5)r(6 - a-3e- z^- z^) 

r(8 - ai3 - 2a245678 - 4e - 2zi - z^- z^^- 2z^ - 2z.) 
r(8 - ai3 - 20245678 - 4e - 2zi - Z2 - Z3 - 2z4 - 2z^ - 2zq) 
r(2 - 0456 - e - 2:4 - Z5)T{2 - a467 - e - ^2 - ^3 - ^4 - ^5) 

r(a45678 -2 + 6 + ^2 + 2:3 + 2:4 + ^5) 
Xr(a4567 + e-2 + Z2 + Z3 + Z4 + 2;5)r(a45678 - 2 + e + Zi + 2:2 + ^3 + ^4 + ^5) 



Xr(4 — ai245678 — 2e — Zi — Z2 — Z4 — Z^ — Zq) 
Xr(4 — 02345678 ^ '^^ ^ ^i — Z^ — Z4 — Z^ — Zq) 

xr{a - 4 + 2e + zi + Z4 + Z5 + Zq) , 



(9) 



We have used this representation to calculate the master planar double box of 
the second type BpL,2{^, ... ,1,0; s,t, rn?; e). The resolution of singularities in e was 
performed similar to the previous cases. The number of resulting MB integrals where 
an expansion in e can be performed in the integrand is again equal to six. We have 
arrived at the following result: 



X 



-BpL,2(l,---,l,0;s,t,m^;e) 
cl{x,y) , c\{x,y) 



^7r^/2e-7Ee 



xy 



:-t) 



l+2e 



+ 



+ Coi{x,y) + Cf^2i^,y) +CQ^{s,t,m ) + 0(e) 



(10) 



where 



cl{x) 



Cl[X) 



In 



^ i^^^y 



-2 In 



+2 In 



1 + a; I + y 

i-y 



l + y 
1 — X 

l + x 



1 1 - a; , 
+ m m 



Li2 

Li2 

1- 



1 — X 

2 
1-y 



- Li2 
-Li2 



X 



l + y 



y 



4 In - + In 

X 1 



1 + a; 
2 

2 

1 — X 



+ Li2 (x) - Li2 (-x) 

+ Li2 {y) - Li2 {-y) 
1 +x 



In 



y 



l + y 



;ii) 



(12) 



Furthermore, 

cl^{x,y) = 4((-mj,+pj,)(3Li3((l-a;)/2) + Li3(l-x) + Li3(-a;/(l-x)) 

-Lia (-2x/(l - x)) - Lis (1/(1 + x)) - Lig {x/{l + x)) + Lig (2x/(l + x)) 

-3Li3 ((1 + x)/2)) + (m, - p,)(Li3 ((1 - y)/2) - Lig (1 - y) - 2Li3 {-y) + 2U^ {y) 

-Li3 (-t//(l - y)) + Li3 {-2y/{l - y)) + U^ (1/(1 + y)) + U^ (t//(l + y)) 

-Li3 {2y/{l + y)) - U^ ((1 + y)/2))) + 41n(-s/m2)((m, - p,)(Li2 ((1 - x)/2) 

-Li2 ((1 + x)/2)) + (m, - pJ(Li2 (-y) - Li2 {y))) 

+2{{py - my) ((2/2 + 4ly -TUy- p y)lA2 {x) + 2(^2 + ^x - "^x)Li2 (-2a;/(l - x)) 

-2{h + h-p^)U2{{2x)/{l+x)) 

-(2 + 2/2 + 4/^ + 4/j, - mj, - 4p^ - pj,)Li2 ((1 + x)/2)) 

+ {m^ -Px){2 - 2/2 -4/j, + m^ - 2mj^ +p^ + 2pj^)Li2 ((1 -y)/2) 

-2Li2 (x) Li2 ((1 - y)/2) + 2Li2 ((1 + x)/2) U^ ((1 - y)/2) 

Hi^x - Px)(2/2 - 4/, + m, + p,) + 2Li2 {x) + 2U2 ((1 + a;)/2))Li2 {-y) 

-((m^ - Px){2l2 - 4/^ + m^ + p^) + 2Li2 (x) + 2U2 ((1 + a:)/2))Li2 (y) 



-2(m, - p^){{h + ly- my)U2 (-2y/(l - y)) - {h + ly- Py)U2 {{2y)/{l + y))) 

+Li2 {-x) {2hmy + AlyTUy " ^J " 2/2^^, - ^lyPy +pI + 2Li2 ((1 - y)/2) - 2U2 (-y) 

+2U2 (y) - 2U2 ((1 + y)/2)) - ((m, - p,)(-2 + 2/3 - 4/, + Aly + m,+p,- 4py) 

-2U2 (x) + 2Li2 ((1 + x)/2))Li2 ((1 + y)/2) 

+Li2((l -x)/2) {{niy -Py){2 - 6/2 - 8/^ - Aly + Gm^ + ruy + 2p^ + py) 

-2U2 ((1 - y)/2) - 2Li2 i-y) + 2U2 (y) + 2U2 ((1 + y)/2))) 

-(m^ -Pa){my -py)\Y?{-s/m^) + {m^my^Wkilx + ly - 1) +7r^))/2 

+ (2(4/3; -m^; -Px){m^ -Px){my - py) \n{- s / m^)) 

+ (8m^(m^ - p^))/3 - GkraKniy - Py) + {imKniy - Py))/3 - 8/^(m^ - p^){my - Py) 



+8l;{m. 



Px)imy - Py) + 2{l2 + m,^)pl{my - Py) - 2pl 



,m,, 



Py) 



-4:ll{m^ + my-p^- Py) + (2(/2 - l)Tx'^{m^ + my - p^ - Py))/3 

+ {6m^my - 2l2{m^ - p^) - p^i^niy + Px))pl - {l^{mx - Px)pI)/3 

-mlirrix - Px) (2/2 + 8/j^ + m^+p^ + 2py) 

+ {Tc'^{myP^ + (m^ - 5p^)py))/6 + ml{2myp^ - 2p^Py + p^) 

+ {2ly{m^ - Px){tt'^ + Qmyp^ + Qm^{my - py) - GniyPy - QpxPy + 12p^))/3 

-All{Px - lyP^ + my{l^ + Px - 1) + Pj/ - IxPy + 2pxPy ^ TU ^{-l +ly- ATUy ^ Py)) 



+ {2lx{my -py){'n'^ 



24ly{mx 



Px 



Qniyp-x - 6pxPy + Qmx{my - Px + Py)))/3 



-SkilyniyPx - (1 + niy)pxPy + nix{{lx - Px)Py + my{px + Py))) 
+4:{-{myPxPy) + mx{~{pxPy) + my{px + Py))) . 



(13) 



The function Cq2(x,|/) is equal to {p^ — mx)bo2{y,x) where the function bo2{x,y) is 
given by the braces in Eq. (9) of jTJj which defines the contribution bo2{x,y) to the 
finite part of the master double box of the first type. 
Finally, 



Cmis,t,m^) 



1 /■! 



—t Jo Jo 



dxidx2 



VI - xi 



Arsh 



X- 



-ty/Xl^/l - X2)l{2m^Xx + X2 - X1X2) 



(4m2 — sx\)x2J [/^vn? — t)xi(\ — X2) + ^vfi^X2 



X (ln(— s/m^) + 2 Inxi — In 4(1 — X\)x2 — sx\[x\ + 0:2 — xiX2)/mf 

1 /■! /■! ln(xi + X2- a;ia;2) 

+ - / / d3:ida;2 , 5 -, 7. r 

s JO Jo \/l — Xi[A.m'^ — sxi)\/l ~ X2[^rn'^ — tx2) 

X (^4 In 2 + 2 ln(-4mVs + Xi) + 2 ln(l - Xi) 
—2 Inxi + 21n(l — X2) — ln(xi + X2 — 3:1X2)) • 



(14) 



We have controlled this result similarly to the previous case, by numerical evalu- 
ation of our finite MB integrals and numerical evaluation by the method of 
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It is not clear whether the two-parametric integrals present in the finite part can 
be expressed in terms of HPL or 2dHPL depending on special combinations of s, t and 
m'^. If the answer to this question is negative one might think about the introduction 
of a new class of functions, e.g., a sort of generalized two-dimensional HPL. In this 
context, let us point out that generalized HPL of various types were introduced in 
These new special functions were defined similarly to HPL, with other basic 



dt 



functions, in particular 1/ Jt(t + 4). 
However, for example, the GHPL 



H{—r, —1; x) 



t{t + 4) 



defined in [T3j equals 



where 



7r\ 1 a TT^ 

2Li2 U, - + - n -2 

^ * ' 3/ 2 18 



y/A + X — a/x 



a/4 + X — y/x 

and Li2 {r,9) = Re Li2 (re'^j is the dilogarithm of a complex argument. Still it is 
not yet clear whether any GHPL can be expressed in terms of polylogarithms and 
HPL. 

From the mathematical point of view it is natural to try to express any new results 
in terms of known functions, which may depend on the initial variables through some 
special combinations. Pragmatically, it is sufficient to present new results in terms 
of some functions which are, probably, new and which can be evaluated with high 
precision at physical values of the given variables. Anyway, the two-parametrical 
integrals presented above can be evaluated numerically in a straightforward way at 
any non-singular point of the variables s, t, m^. 

2.3 Non-planar double boxes 

For the general non-planar double box (see Fig. Ic), one can derive the following 
eightfold MB representation 



BNp{ai,...,as;s,t,u,m ;e 



.in^/^)^(-lY 



n,=2,4,5,6,7r(a,)r(4 - a4567 " 2e){-s)-'^+^^ 

r(a5 + Z2 + Z4^)T{ar + Z3 + Z4)r(4 - 0455677 - 2e - Z2 - Z3 - 2zi) 
r(ai - Z2)r(a3 - 23)r(a8 - Zi) 



r(2 — 0567 — ^ — Z2 — Z4 — z^)T{2 — 0457 — e — Z3 — Z4 — z^) 

r(4 — 0455677 — 2e — Z2 — Z3 — 2^4 — 2^5) 

T{a8 + zi- z^ + Z7)T{8 - ays - 2a245678 - 4e - 2:2 - 23 - 2^5 - 2^7 - 2^8) 



X 



X 



r(6 — a — 3e — 2:5) 

r( — gg — Zi + ^4 — 2:7 — Zg)T{A — 02345678 — "26 — Z2 — Z^ — Zq — Z^ — Zg) 

r(8 - ai3 - 2a245678 -4:6- Z2- Z3- 2Z5 - 2zq - 2Z7 - 2Z8) 

r(4 - ai245678 -2e- Z3- Z5- Zq- Zy - Zs)T{a28 + Zi~ Za + Zt + Zs) 

r(a245678 -2 + e + Zi+ Z2 + Z-i + Z^ + Z7 + 2Z8) 

xr(a4567 + e - 2 + 2:2 + 2:3 + Z4 + ^5 + Z8)V{a - 4 + 2e + 2:5 + 2:6 + ^7 + Zs) 
Xr(a245678 -2 + e + Zi+ Z2 + z^ + z^ + 2zt + 2z8) . (15) 

Although the number of integrations is rather high one can proceed also in this case. 
However, it turns out that the massive non-planar case is rather complicated. We 
have performed the resolution of singularities for the non-planar master planar double 
box. Here is our result for its double-pole part in e: 



d/2\ 

stu 1 + X 1 + z e'^ ' Ve 



I trc ' ] xz \ — X 1 — z I f\ 

Bnp{1, ...,l,0;s,t, m'; e) = ^ ^ In ^— In ^^ 7^ + O ( - 



where z = l/Jl — ^lW? /u. A one-parametric integral which can hardly be expressed 
in terms of known special functions is present already in the 1/e part. 

Unfortunately, we are unable to check our preliminary result (for the coefficient 
of the simple pole and for the finite part) numerically by the method of J2SI for 
the following reason. As in the case of the massless non-planar double boxes (see 
the second reference in j2]), it is natural to treat the Mandelstam variables as not 
restricted by the physical condition s + t + u = 4m^, because this condition does not 
simplify the calculation. Then the natural procedure is to perform the calculation 
for some extension of the given Feynman integral as a function of the variables s, 
t, u (and, now, rn?) to the Euclidean domain, with s,t,M < 0, and perform the 
analytical continuation to the physical domain in the result. When one starts from 
the alpha/Feynman parameters of the whole diagram one naturally arrives at an 
extension of the physical condition just by considering u as an independent variable. 
Since, in the massless case, a fourfold MB representation was derived starting from the 
global Feynman parameter representation, the same extension as in i23^ was implied 
there. In the massive case however, if we start from Feynman parameters and try 
to separate various terms entering functions present in the integrals over Feynman 
parameters, we do not see any possibility to arrive at a MB representation with 
the number of integrations less than ten. The eightfold MB representation ()15|) was 
derived in another way, by introducing some Feynman parameters upon integrating 
over the first loop momentum, then integrating over the second loop momentum and 
completing the procedure of introducing MB integrations. After this procedure, the 
variable u is considered as an independent variable. It turns out that this extension to 
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non-physical Euclidean points defined by p5|) differs from the version implied within 
P^ . based on the global Feynman parametric representation for the initial diagram. 
This can be seen by checking (I15|l in limiting cases where some of the indices Oj 
are zero and where one has explicit analytical results: In some of these limits, it is 
necessary to use the physical condition s + t + u = Am? to prove agreement with the 
known result, while the agreement is achieved without this condition when one starts 
from the global Feynman parametric representation. 

Thus, to perform numerical control of our results, a generalization of the method 
of sector decompositions 23j to points with kinematical invariants of different signs 
is necessary and will be developed in the near future. Then one will be able to make 
reliable results for massive non-planar double boxes. 
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